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SUMMARY 


Instead of solving the nonlinear differential equation vhich governs 
the compressible flov, an approximate method of solution by means of the 
variational method is used in this report. The general problem of steady 
irrotational flow past an arbitrary body is formulated. Two examples 
were carried out, namely, the flow past a circular cylinder and the flow 
past a thin curved surface. The variational method yields results of 
velocity and pressure distributions which compare excellently with those 
found by existing methods. These results indicate that the variational 
method will yield good approximate solution for flow past both thick and 
thin bodies, at both high and low Mach numbers. 


INTRODUCTION 


The advance in aerodynamics has been greatly accelerated by the 
interpretation of good experimental results in the light of theoretical 
predictions. When the flow velocity is small, the change in density of 
the fluid in the flow is so small that the fluid may be considered to be 
incompressible. Many useful theoretical results have been obtained from 
a consideration of the two-dimensional irrotational flow of an incom- 
pressible perfect fluid. When the flow velocity becomes large, the 
change in density of the fluid in the flow can no longer be considered 
small and the flow must be considered to be compressible. While the two- 
dimensional incompressible flow can be described by the well-known 
Laplace equation, .the corresponding compressible flow. can only be 
described by a complicated nonlinear differential equation. Except' for 
a few simple cases, the exact solution of this differential equation is 
still unknown. . . 

In the case of two-dimensional irrotational subsonic flow, there 
are at least three different methods of small perturbation or lineariza- 
tion. The small -perturbation methods have been carried out by developing 
the velocity potential in terms of the Mach number or of a thickness 
parameter. Obviously, the development in terms of the Mach number limits 
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the applicability of the method to flows at small velocities and the 
development in terms of the thickness parameter limits the applicability 
of the method only to flows about very thin bodies. Methods for calcu- 
lating the second- and the third-order approximations in these methods 
have been carried out in some cases/ but the increased mathematical 
difficulty with each approximation prevents their application to the 
higher orders of approximation. Another method of linearization is the 
so-called hodograph method wherein the nonlinear differential equation 
is transformed into an exact linear differential equation, thus removing 
the difficulty of the nonlinearity of the governing equation. However, 
in doing so, the fulfillment of the boundary conditions is made much 
more difficult. The application of this method to actual problems 
requires some simplifying ass\imptlons which can be shown as being only 
approximately correct . ' 

The variational method, as carried out in this report, on the other 
hand, can be used to study the compressible flow without linearization 
and the computation can be made to obtain a solution of high accuracy 
without any mathematical difficulty. The variational method has also 
the eidvantage of being a direct method by means of which one may be able 
to calculate the flow passing a given airfoil at different Mach numbers. ' 
In the variational method, instead of solving the nonlinear differential 
equations of motion, a variational integral is first formulated and 
approximate solutions are then calculated by the Rayleigh-Ritz procedure . 

Two numerical examples are carried out in this report to ascertain 
the applicability of this method. They are for the two-dimensional 
irrotational flows past a circular cylinder and a thin bump as proposed 
by Kaplan. (See reference 1 .) These examples are chosen because in the 
case of circular cylinders similar results by the Rayleigh-Janzen method 
are available, and in the case of the Kaplan bump results by means of 
the perturbation method in terms of thickness parameter are available. 

A comparison of the results by the present method with those computed 
by the other methods shows excellent agreement. Since the Rayleigh- 
Janzen method is known to give good results for thick bodies at a rela- 
tively low Mach number and the perturbation method in terms of thickness 
parameter yields accurate results for thin bodies up to high Mach numbers, 
the results of this report indicate that the variational method gives 
good approximate solutions for both thick and thin bodies and at both 
high and low Mach numbers. 

This Investigation was carried out at the Daniel Guggenheim School 
of Aeronautics, College of Engineering, New York University, under the 
sponsorship and with the financial assistance of the National Advisory 
Committee for Aeronautics. The author wishes to acknowledge his indebt- 
edness to Messrs. J. R. Knudsen and H. H. Hilton who carried out most of 
the numerical, calculation. Grateful acknowledgment is hereby also made 
to Professor K. 0 . Friedrichs for his kind encouragement and inspiring 
discussions and to Professor F. K. Teichmann for his understanding 
cooperation. 
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SYMBOLS 

velocity of soxmd 
pressvure coefficient 

pressirre coefficient for incompressible flow 

variational integral 

Mach number (q[/a) 

inward normal direction 

pressure 

velocity vector 

magnitude of velocity vector 

maximum possible velocity of (flow 

position radius in flow field 
velocity component in x-direction 
velocity of undisturbed stream 
velocity con^pohent in y-direction 
rectangular coordinates of a point in the fluid 
coordinates of curved surface divided by semichord of shape 
+ iy = re^® 

thickness parameter 

constant defining thickness of Kaplan' s bump 

ratio of specific heats at constant pressure and constant 
volume 

angular position 


0 


velocity potential 



k 
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p mass density of fluid 

I ,T] rectangular coordinates in plane of curved surface 

I 

C = I + iT] 

r strength of circulation 

Subscripts: 

x,y differentiation in corresponding direction 

o conditions of undisturbed stream 


THEORY 


For steady two-dimensional irrotational flows, there are fovir hydro 
dynamical variables; The pressure p, the density p, and the two 
components of velocity u and v. All these are to be determined as 
functions of the coordinates x and y by the following four 
equations: 

(1) The equation of continuity (one equation) 

div (pg) = 0 , • 


or 


^pu , 3 pv 



( 1 ) 


where q is the velocity vector. 

( 2 ) The equations of motion (two equations) 
^ ^ -1 St 


du , du -1 dp 


Sv dv -1 dp 

^ SE ^ ^ = -P ^ 

(3) The equation of state (one equation) 

p = A + Bp^ 


Sv 


-1 Si 


( 2 ) 




( 3 ) 


where A/ B, and 7 are constants and the significance of A and B 
will be discussed later. 
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Most aerodynamic problems are concerned with either flow from rest 
or flow which is parallel and uniform at infinity. Both of these types 
axe irrotational to begin with^ and will remain so provided that fric- 
tion and shock waves are neglected. Thus, irrotational flow is of great 
interest in aeronautical applications. 

The condition of irrotationality is 

rot q = 0 
or 

^u Sv _ n 

This condition can be satisfied by introducing a velocity potential 0 
such that 



that is, q is a gradient of 0. 

For the velocity field given by equations (k) to represent a physi- 
cally possible field of flow, it is necessary that the continuity equa- 
tion (1) be also satisfied. By means of equations (l), ( 2 ), anfl (3), p 
and p are eliminated, and one obtains the following fundamental dif- 
ferential equation: 

, - 0y^0yy - 20x0y0xy = 0 - (5) 

in which 

where a is the velocity of sound in the fluid and the subscript o 
refers to quantities at infinity. Equation (5) is a nonlinear differ- 
ential equation, and its exact solution is found to be extremely diffi- 
cult. Instead of solving the differential equation exactly, an approxi- 
mate solution may be obtained by means of the variational' method. In 
this method the first step is to formulate a variational principle. 
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A vaxiational principle for the compressible fluid was first 
formulated by Hargreaves. (See reference 2.) A somewhat more detailed 
discussion of this variational principle was made later by Bateman. 

(See reference 3 .) A simplified proof of the principle was given by 
Wang. (See reference 4.) It seems, however, that the following formu- 
lation is more complete from the mathematical point of view anH is 
therefore included here . 

Consider the following variational integral 



where the integration is taken over the given domain D. Varying p, 
u, V, and 0 independently and setting the first variation equal to 
zero, one obtains 



= 0 


where S is the boundary of the domain D and q^^ is the velocity 
component which is normal to the boundary S. ' 


The Euler's equations of the integral I, according to the calculus 
of vaxiatlons, are then 


f'(p) 



+ |(u^ + v^) = 0 


(7) 



( 8 ) 



(9) 


Spu Spv _ 0 

sr ^ 


(10) 
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and the natural 'bomdary condition must be such that 



pB 0 qn dS = 0 


( 11 ) 


Equation (lO) is the continuity equation (l) and equations (8) and (9) 
are the same as equations (4)^ which implies that the condition of irro- 
tationality is satisfied. If the pressure p is given by the equation 


p = -pf '(p) + f(p) (12) 

equation ( 7 ), when combined with equations ( 8 ), ( 9 ), and ( 10 ), indicates 
that equations ( 2 ) are satisfied. 

Thus the variational principle ( 6 ) leads to the fundamental equa- 
tions of a two-dimensional compressible fluid if p is of -the form 
given by equation ( 12 ) and the condition (ll) is satisfied. Equa- 
tion (i2) is a differential equation of the Clairaut type and its solu- 
tion, when p = A + Bp7, is 

f(p) = A + q ^^2 p/2 _ Bp7/(7 - 1 ) - (I3) 


tion ( 11 ) is satisfied if the domain is finite because the boundary con- 
^ilions in fluid dynsmics are either that- the boundaiy is a stresm sur- 
face, thus q^^ = 0 , or that the velocity is prescribed; that is, 50 = 0 . 
For the domain which extends to infinity, since a part of the boundary S 
is infinite, the integral may result in a finite value, although the 
integrand may approach zero at infinity. This case has already been dis- 
cussed in detail by Wang (reference 4) and in such cases this finite 
value of integral (ll) should be subtracted from integral ( 6 ). 

In applying the Rayleigh-Ritz method to the problem, it is more 
convenient to use only one variable in the variational integral, say the 
velocity potential 0 , and to express the other variables u, v, ' 
and p in terms of 0 by means of equations ( 8 ), ( 9 ), and (lO). In so 
doing, equations ( 8 ), ( 9 ), and (lO) are automatically satisfied and the 
first variation of equation ( 6 ) then leads to the fundamental differential 
equation ( 5 ). With equations ( 8 ), ( 9 ), (10), and (12) satisfied, it can 
be easily verified that the integrand of I is just the pressure p. 

The Rayleigh-Ritz method consists essentially of choosing a serios 
of functions which satisfies the boundary conditions but has undetermined 
parameters, such as A^j, i,j = 1, 2, . . ., in the functions. If the 
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variational integral is denoted by I, then by substituting these 
functions into I and setting the first VEiriation of I equal to zero, 
the following set of simultaneous algebraic equations is obtained: 






= 0 


■ij 


= 1 , 2 , 


The parameters A^j can thus be determined by solving these simultaneous 

equations. The Rayleigh-Ritz procedure, therefore, gives a solution 
which satisfies the boundary conditions exactly but satisfies the dif- 
ferential equations only approximately. 

Consider the case of a compressible-fluid flow past a circular 
cylinder with unit radius. The velocity potential 0 in this case may 
be assumed tobe- 

0 0 q + 02 . ( 1 ^) 

with - 

01 = u(r + j)cos 0 + 0 (15) 

and • ■ , 

‘ £ £ fc ■ - V “«) . (16) 

where T is the strength of circulation, A ,„„ and are undetermined 

parameters, and r and 0 are the polar coordinates. It is evident 
that for the assumed value of 0 the boundary conditions are satisfied; 
that is, at r = w and q = U, U being the velocity of the undisturbed 
flow at Infinity and on the surface of the body, 80/8h = 0 or the 
normal component of the velocity is zero, n being the inward normal 
direction. 

The circular cylinder in the complex z -plane may be mapped into 
cylindrical bodies of arbitrary shape in the complex ^-plane by the 
mapping function ^ = f(z). When the mapping is conformal, it was proved 
in reference 4 that 0(0 still satisfies the boundary conditions. 

With 0 as assumed in equations ( 15 ) and ( 16 ), the final variational 
principle (reference 4) may be written as follows: 
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I = 


rr 

2 c 

l(z)^ 

7/(7-!) 


ij 

Q- 

max 

ld ^2 


dz 



dzl_ 




r dr d© + 


z -plane 
2 r 


- 1 ISm! 


nax 


-u^) 




UjtA 


■11 


(17) 


where the last term is due to the fact that the integral (11) is not 
zero and 



The final velocity in the plane of this arbitrary body is 


= q(z)^ 



(19) 


The integral equation (I7) can render itself into analytical inte- 
gration only when 7/(7 - 1 ) is an integer. For different values of 
7/(7 - 1), the corresponding values of 7 are as follows: 


7/(7 - 1 ) y 


2 2 

3 1.5 

4 1.333 

The isentropic value of 7 for air is 1 .^ 05 , which is between the 
values of I.5 and 1.333. Among these values of j 7 = 1.333 gives 
the closest approximation to the isentropic value. However, in this' 
case the integrand contains an expression which is raised to the fourth 
power and the labor in carrying out the integration becomes excessive. 
In order to simplify the numerical work, 7=2 is used in the subse- 
quent calculation. Chaplygin and Von Karm^ and Tsien used 7 = -1 in 
their solution by the hodograph method. The justification of taking 
7 = -1 has also been discussed by the latter two-authors. The idea is 
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that the velocity in the flow field does not differ too nnach from the 
undisturbed velocity and therefore only a small portion of the p - p 
curve will be used. If one makes the curve with 7 = -1 tangent to the 
curve with 7 = 1.405 at the point that represents the. conditions at 
infinity (pq,Pq) and uses the curve with 7 = -1, the error will not be 
too large. The use of 7=2 Instead of 7 = -1 gives even better 
approximation as can be readily seen from figure 1. Although the curves 
with 7 = -1 , and 7 = 1.405 may be made tangent to each other at 
(Pp^Po)^ second derivatives have opposite signs. When the curves 

with 7=2 and 7 = 1.405 are made tangent to each other at (po^Po)^ 
however, a small portion of the curves is very close together, and thus 
a much better approximation is provided. 

The isentropic pres sure -density relationship is as follows: 



( 20 ) 


According to the previous discussion, the constants A and B in 
equation. ( 3 ) can be determined from the condition 'that the curve given 
by equation ( 20 ) and that by equation ( 3 ) are tangent to each other at 
(Pq^^^o)* This yields two conditions) namely. 


and 



( 21 ) 




( 22 ) 


Thus in the case of 7 = 2 , the constants A and B are found to be 
A = 0.2975 Pq and B = 0*7025PQyp^ and equation (20) becomes 

0.7025p„p^ 

P = 0.2975Pq + ( 23 ) 

Pq 

For 7 = - 1 , the corresponding pres sure -density relationship results in 


p = 2.405p^ 


1.1.05p^p^ 


P 


(24) 
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For flow past the circular cylinder, the range of variation of p/p - 
is approximately from 0.2 to 1.2. In this region, it can be easily^ 
seen from figure 1 that equation (23) does indeed give a very good 
approximation of equation (20). 


FLOW PAST A CIRCULAR CYLINDER 


The problem of determining the flow past a circular cylinder with- . 
out circulation has been carried out by Wang and reported in reference 4. 
A detailed discussion has been made there for the flow at Mq = 0.4. 
Because of the symmetry of the problem, 0 may be assumed as follows: 

0 = U (r + ijcos 0 + A^^^i - ^jcos 0 + - ^jcos 30 + 


In expression (25), six unknown parameters are taken. The variational 
integral is 



(26) 

where 7 is taken as 2. 


Substituting 0 as given by equation (25) into equation (26) and 
differentiating the resulting integral with respect to the parameters, 
one obtains 


ai/SA^^ = 0 

= 0 

si/aAji 

= 0 

= 0 

^I/^l 
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Integrating these equations, there results 

(2-07k07q^^ - 7.31852u2)aj^ .* 7 . 0666 tu^A ^2 + (0.325926(3^^^ - 
1.03TOirtj2)A3i + 0.915555 u2a 33 + (o.ll4286qj„g^2 . 0 . 3^*Ol36u2) 

7.0666 Tu 2 a 3 _i + (l3.3333qniax^ - + 1 . 20000 u2a33_ + 

(l.60000(]jjiax^ - 5.9733 ^u^)a 33 + 32 . ii-W-lrtJ^Ai^ + 0 . 4iwD8l6u2A5]_ = K2 

(0.325926(3^^2 _ i.03704u2)a^^ + 1.20000 u 2 a ^2 + (o . 09007^1 q^^2 _ 
0.207U07U^)a31 + o.192254u^A 33 + (o.oiai428qjjjg^^ - o.o8io885U^)a5i 

0.915555u2a3_i + (l.60000Qj^g^2 _ 5.97334u2)ai3 + 0.192251|u2a33_ + 

(o.277333qmax^ - 1 . 01l8lU^ ) a 33 + k.92kkku\^ + 0.0813062U%1 = K4 

i 

, 32.mku\^ + 1^.92444u^A 33 + (35.8519qmax^ - 120.948 o^)ai5 = K 5 
(0. 114286(3^^^ - 0.340136 u^)a^^ + 0.44o8l6u^A^^ + (o.o4ll428(3^^^ - 

0.0810885U^)a 23^ + 0.0813062U^A^2 + (o.02l8309q ^^^ ^ - 
0.0360739U^)a^^ = Kg 

' ' } 

where K]_, K 2 , K3, Ki^., K^, and Kg are s-ums of nonlinear terms in 

■^11-' ^13^ ■^31^ ^^33^ ■^15'> -^1 their complete expressions 

are given in reference 4, except that. the first' term 'in Kq should 
be 4.22222U3 first term in K3 should be -6.66666U3; the 

corresponding values given in reference 4 axe due to misprints. These 
equations can be solved as follows: First, assume the values of the A's 

to be zero in the K's. Equation (27) then becomes a system of linear 
simultaneous equations which can be solved by the method of solving 


= K- 


>(27) 
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linear simultaneous equations suggested by Grout (reference 5 )* The 
values of A-| 1 , -^13^ -^31^ ^33^ ^15^ ^1 computed may be 

taken as the first approximation. Use the values of the A's so deter- 
mined to compute the K's and solve for these parameters again. This 
gives the second approximation. Repeat the cycle until the desired 
accuracy is obtained. To make the convergence more rapid, after the 
third approximation, instead of following the above-outlined method 
rigidly, one may extrapolate the results and use the values of the A's 
so determined to repeat the process. This alternative method is found 
to be especially time saving in the computation for the cases where the 
Mach number is relatively high. A detailed example of the computation 
may be found in reference k. 


The parameters for flows at various undisturbed-stream 

Mach numbers are computed and are tabulated in table I. In order to 
see the convergence of the series, these parameters are computed by 
taking one parameter only, two parameters, three parameters, and so 
forth, up to six parameters. The maximum velocities at various Mach 
numbers are tabulated in table II for all the six cases. The convergence 
of the series is seen to be satisfactory except at Mq = 0.5 where the ' 
series is apparently divergent, and with four terms the maximum velocity 

is already greater than the maximum allowable velocity . For 7=2, 

P Tnax 

/qjnaxV 2 

— J = 2 + Mq , The velocity distribution over the cylinder at 


Mq - 0.4 conq)uted by the variational method is plotted in figure 2 as 
are the results by the Rayleigh-Janzen method to the third approximation 
as computed by Imai (see reference 6). At. various Mach numbers of the 
undisturbed flow, the maximum local Mach number is plotted in figure 3 
and the local Mach mmiber over the surface of the cylinder is plotted 
in figure 4 . The velocities computed by the Rayleigh-Janzen and the 
present method show excellent agreement. 


In order to study the effect on the velocity of using 7=2, the 
next value closer to the isentropic one, that is, 7 = I. 5 , is now taken 
to compute the velocity. Taking only the first term in the velocity- 
potential series, the equation Si/SA3_q^ = 0 is obtained as follows: 

- 3.28889u 5) - (o.5l85l8q^^^ - 3 . 65926 q^^V ^ 

5.928o4U^)a3_i - (3.21005U^ - 0.82222qjjj^^juAi;L^ - ^1.05T66U^ - 

o.ioi59q^^^)A^^3 - 0 .i6490UA^^^ - 0.010564a^^5 = 0 (28) 
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The maximum velocity with 7 = I.5 is also^included in table II for 
cou5)arison. It is seen that these velocities are smaller than the cor- 
responding ones computed by taking 7=2. However, the velocities also 
increase with increased number of terms in the 0 series. Since o nl y a 
finite niimber of terms can be taken in the numerical calculation, the 
velocities computed will always be smaller than the corresponding exact 
values. In such cases, the effect of taking 7=2 is to increase the 
velocity, which compensates to some extent the effect of taking only a 
finite number of terms. 


The pressure coefficient at a point whose pressure is p may be 
defined as follows: 



For 7=2, 



The pressure coefficient is computed from the 'above formula and is 
plotted in figure 5- Results by the Rayleigh-Janzen method and the 
variational method compare very well, with the former yielding slightly 
lower results at small angles. The pressure coefficient by the Rayleigh- 
Janzen method was computed by taking 7 = 1.405 in equation (29), which 
becomes upon substitution: 


^ 0.7025Mq2 

Both the variational and Rayleigh-Janzen pre'ssure-coefficient values are 
lower than those obtained by the Von Karmdn-Tsieri method and generally 
higher than the Prandtl-Glauert approximation. For further comparison, 
the pressure -coefficient ratio is plotted in figure 6 where 

Sa pressure coefficient of the corresponding incompressible 

flow. 


1 + O.2025M0" 



13.469 1 

(‘■^1 



(31) 
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FLOW PAST A CURVED SURFACE 


In order to obtain further evidence of the applicability of the 
method, the compressible flow past a thin curved surface will next be 
considered. This example is carried out for two purposes. First, from 
the previous example , the variatio.nal method is seen to give very good 
results for flow past a thick body where the Mach number is relatively 
lowj it is to be ascertained now whether the method will also give good 
results for flow past a thin body at high Mach numbers. Second, there 
may be _some difficulty , in carrying out the method when an arbitrary 
mapping function is used. If this is so, a study should be made to 
remove such difficulty. 

The compressible flow past a thin bump has been carried out by 
Kaplan (reference 1) by the perturbation method in terms of a thickness 
parameter. Results are given for flow past a series of bumps of various 
thicknesses. Kaplan's bmp may be obtained by mapping conformally a 
circle of unit radius in the complex z-plane into the complex ? -plane by 
the mapping function 



(32) 


Again take six parameters in ^ as in the previous example, noting 
that z = re^®. The equations = 0 may now be written in the 

following form: 



- =0 (33, 



l6 
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where 



obtained from the mapping function (32). The quantity = 1 when 

i = j = 1 and 0 for other values of i and J. 

Substituting into the integrals (33) the expressions (25) and 
carrying out the integration first with respect to the angle 9, the 
resulting integrals are of the type ’ , 

■ ' P 2jt 

Xf(sin 6 cos 0)d0 

0 


This integration can be somewhat simplified by the use of the theory of 

residue. Let Z = e and express the integrand of I as a function 
of Z. Each I then consists of a sum of integrals of the two 
following types: 



(1 + d^)[r^ + - l) 


and 
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There are altogether 66 integrals encountered in the calciila- 

tion, and they axe given in appendix A. With these expressions of In 
substituted into the integrals, the resulting integrands axe rational 
functions of r with coefficients in d. It is observed that the final 
coefficients of the unknown parameters Aj^j axe sums of the following 

functions : 


Y 


n 




r 


n 


dr 





where n may be either positive or negative. There are altogether 20 
such Yn expressions encountered in the calculation, and they axe 

given in appendix B. It may be noted that terms like [^ogg (l - 

and (logp r) _ appear after the integration. These terms however 

will be canceled by each other in the final summation. This fact may 
conveniently be regarded as a check to the correctness of the computati 


2 

The numerical values of d may be substituted at this stage befo; 
the final integration with respect to r is caxried out; otherwise the 
computation will be extremely laborious. The integrands now may be tafc 
as the products of r expressions and In and then the final euiswers 
can be immediately written down from the known Yn values. There axe 
eight r expressions involved, and they are 


P = r + - 
r 


^ r 3 
3r^ 


R = 


3r^ 5r^ 


S = 


7r' 


K = 1 - ^ 
r^ 


w 1 1 

M = r + —r 
r^ r^ 

H 


r** 

1 


P 

As a numerical example^ d = 0.075 is taken^ which corresponds to 
a bump with a maximum thickness of 5*13 percent of the chord. The shape 
of the bump is shown in figure 7* d^ = 0 . 075 ^ values of Yn 

axe calculated and included in appendix B. After dropping a common 
factor of 4, integrals (33)^ which axe equivalent to 


si/aAii = 0 

31 / 3*33 - 0 


Bi/SA^5 = 0 


= 0 
= 0 
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become 

(o.5l85l9q^^^ - 1.193255 U^)a^^ + 0 . 55^586u^A^^ + ^0.0814815(3^^^^^^ - 
0.i69149U^)a2^ + 0.0285596u^A^ 2 + 0 . 2l6663U^A^^ + ^0.02857l4g ^^^ ^ - 
0. 05468890 ^) a^jl = % 


o.554586u^A3^^ + ( 3.333333(3 ^^^ - 8.i88034u^)a^2 + 0.112085U^A^^ + 
(o.4ooooo q^^/ - i.012468u^)a33 + 2.347577u\ 5 + o.o457i58u^A^^ = 


(o.08l48l5(3jnax^ - 0.169149 U^)aii + 0.112085u\3 + (o . 0225185(3 j^qx^ 
0.0420399U^)a3i + 0.00795193U^A33 + 0.0706980U^Ai5 + 
(o.0102858qj„a^^ - 0.0178331 U^)a5i = K 3 


>(34) 


0.0285596U^Aii + (o.400000qmax^ - 1. 0124680^ jAi3 + 0 . 00795193U^A3i + 
(o»o693333cUfly^ - 0 .i676i6u^)a 33 + 0 . 3198l8U^A3_5 + 0.00397346u^A5i = 


o. 2 i 6663U^A3^3^ + 2. 3475770 ^ Aj ^3 + o. 07 o 698 ou^A 3 i + o. 3198 i 8 u^A 33 + 

(8.962963q„^^ - 22.05199 U^)ai 5 + 0.0356363U%i = K5 


(o.02857l4qjjjg^^ - 0. 05468890 ^ ^A^^ + 0.0457158u^a^3 + 
(o.o102857(3 ^„/ - ,0.0178331U^)a 3^ + 0 . 00397346u^A33 + 
0.0356363U^A^^ + (o-0054577(3^^^ - 0.00887571U^)a^ 3_ = Kg 

The Kjj expressions are found, in appendix C. 
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These equations ( 3 ^) are solved again by the method of successive 
approximation as outlined in the previous section. At relatively high 
Mach numbers, however, the method as outlined before becomes oscilla- 
tively divergent. In such cases the average values of Aij between two 
consecutive cycles of computations may be used for the next approxima- 
tion. The process is then again convergent. Such a method was also 
found necessary by Wang in a different problem where another set of 
cubic simultaneous equations was involved (reference 7 )* 

At Mq = 0.9, the method described above again failed when six 
parameters were used. Many schemes of computation were tried and it 
was not possible to calculate the values of the parameters by methods 
of successive approximation. It was possible, however, to calculate 
these values when up to four parsuaeters were taken. This indicates that 
at this Mach number Irrotatlonal solution may not exist (reference 8 ) 
and shock waves have possibly occurred. It can also be seen that the 
velocity computed by Kaplan begins to diverge at about Mq = 0 . 9 . The 
velocity series computed by Kaplan becomes in this case 


Q 2 

y = 1 + ^ 3 ^ + • . • 

= 1 + 0.17614-72 + 0.050675 + 0.059108 + . . . 


where t is the maximum thickness ratio and is equal to O.O51282 in 
this case. 


The values of Aj^j^ag are given in table III at Mq = 0 . 5 ^ 0 . 75 > 

0.83, and 0.90. The velocity distribution over the surface of the 
b\imp is given in table IV and plotted in figure 7 to compare that found 
by the variational method with the results by Kaplan. The agreement is 
very good. At Mq = 0 . 9 ^ the maximum velocities q/U with' one, two, 
three, and four parameters are 1.1179> 1.2628, I.2716, and I.263O, 
respectively. The successive values of the velocity do not show a sign 
of divergence in this case. The velocities at various points on the 
surface of the bump are obtained as follows. For the mapping func- 
tion (32), the shape of the bump in the C-plaiie is given by the following 
parametric equations. 


1=2 cos 9 - ^(3 cos '9 - cos 36 ) 
^2 

Tj = ^(3 sin 0 - sin 30) 


s 






(35) 
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and the maxi- 


The length of chord is given by (2|) or (2 ~ 

®=0 o' 3 / “ 

mum thickness is given by or' kdV3. The thickness ratio t 

is then equal to 2d^/(3 - d^). If the unit length is taken as the 
semichord of the shape, equation (35) becomes 


X = cos 6 - j^(cos 6 - cos 30) 


Y = ^(3 sin 6 - sin 36) 






For different values of X, 9 may be found from equations (36). 
the surface of the bump. 


(36) 


On 


III = ||l - (1 - d2)cos 20 - d^ cos - d2)sin 

- 


20 + 


d sin 40 


1 

12 


(37) 


The velocity is then given by the following expression: 
_/l 30 ^ 


3 , = 
U 


^®/r=l 


IT 

dz 


dz 


^13 . ^*^^33 
^oMo , 


^^0^0 


^ 0 . 6666'IAi2_ ^ 0*133333A32. ,0.057142A52\ 

^0^0 VV a^Mo 7 

3-33333Ai3 


Bin 0 + 


I sin 30 + 


^0% 


sin 50 


(38) 


For different values of X, the corresponding values of sin 0, sin 30j 


sin ^ 9 , and 


di I 

^ are given in table V. 


The pressure coefficients at various Mach numbers axe given in 
table VI and are plotted in figure 8. The pressure coefficients com- 
variational, the perturbation; the Von KfitnuEfii-Tsien, and 
the Prandtl-Glauert methods are given in table VII and are plotted in 
figure 9 for comparison. 
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RESULTS AMD DISCUSSION 


Compressible flows past a circular cylinder and a thin curved sur- 
face have been computed by the variational method. In the case of flow 
past a circular cylinder, the velocities and pressure coefficients com- 
puted axe compared with the results obtained by Tmai who used the 
Rayleigh-Janzen method up to the third approximation. In the case of 
flow past a thin curved surface, the velocities and pressure coeffi- 
cients axe compared with the results obtained by Kaplan who used the 
perturbation method in terms of a thickness parameter up to the third 
approximation. The agreement of the results between the present method 
and the existing methods is very good. It thus shows that the VEirla- 
tional method will give good approximate solution both for thick and * 
thin bodies and both at low and high Mach numbers. 

In carrying out the variational method, 7 is taken to be 2 instead 
of the Isentropic value which for air is 1.405. A study is made by 
taking 7 = I .5 in the case of flow past a circular cylinder. The 
results show that the, effect of using a larger value of 7 is to 
increase the maximum velocity slightly. The same conclusion may be 
obtained by a study of the results of the Rayleigh-Janzen method. 
Theoretically, if one can take a complete set of functions in the 
velocity-potential series, exact solution may be obtained. In actually 
carrying out the method, however, only a finite number of terms can be 
taken. Although the , introduction of more terms into the solution does 
not involve any additional mathematical difficulties, each new term 
appreciably increases the computational labor. With six terms, it is 
observed that the maximum velocities increase with an increased number 
of terms. Thus the slightly higher values for maximum velocities axe - 
increments in the right direction. By using a modified pressure-density 
relationship instead of the Isentropic one, the pressure coefficients 
axe haxdly affected by ein increase in 7. 

When the undisturbed velocity increases to a value at which shock 
waves will probably occur, the nonlinear terms of the imdetermined 
parameters axe no longer small compared with the linear terms. The 
method of successive approximation in solving these simult^eous equa- 
tions suddenly fails.' None of the many schemes tried worked when more 
than four parameters were taken. An, investigation has been carried' 
out to study the reason why it is so. It was foxmd (reference 8) that 
the irrotational flow solution fails to exist when the Mach number is 
increased beyond a certain limiting value. For the cases computed, 
it appears that this limiting value decreases as the number of terms 
Increases. This explains why the method of successive approximation 
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fails to work with more parameters while it was still possible to 
obtain a solution with fewer parameters. 


New York University- 

New York, N. Y., June I7, I949 
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APPENDIX A 


THE In EXPRESSIONS 


The In expressions encountered in the calculation of flow past a 
curved surface by the variational method are as follows: 

p 2 it 

lx = X cos^e d 0 

J 0 


Ttr 


4(1 + d 2 )(r^ + d 2 )(r^ - d'^)(r^ - l) 

(1 + - 4 d^ + d^)r^ - d^(l + d^Yj 


[3(1 + d^)r® + 4(1 - d^)r^ + 


I2 = / X cos^0 COS' 30 d0 
^0 


itr 


[[1 + d 2 ) 




3(1 - d^)r® + 


4(1 . d2)(r^ . d")(r 4 . . 1) 

(1 + d 2)(3 - 4 d^ + 3 d^)r^ + (1 - d®)r^,- 3 d^l + d^jr^ - d^(l - d^^ 


p 2 rt ■ 

lo = / X cos'^G cos ,5^ d 0 
Jq 


R 

1 >»10 , ofl i 

4(1 + d^)(r^ + d^)(r^ - d^)(r^ - l) ^ 

1- Q 7 + 311 + d ) 


3(1 - d®)r^ + (1 - 3d^ - 3 d^ + dl°)r^ - 3d^(l - d^)r^ d^(l + d^) 
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'~k ~ I ^ cos 6 cos 2iQ de 

J 0 


Ii. = 


P ^ 2(1 + d2)rl2 + 2(1 - d^)rlO + 


4(1 + d 2 )(r'^ + d 2 )(r^ - d^)(r"^ - 1) 

(1 + d2)(l - 3d2 + d^)r8 + 2(1 - dS)r^ + (1 - d^)(l - d^)r^ - ' 
2d^(l - d^)r2 - d^(l + d^) 


>2rt 


Ir = 


Lc = j X COS 6 cos 30 cos 5® d6 
^0 


+ d 2 )rl^ + 2(1 - d^)rl 2 + 


4^1 + - l)r^ 

(1 + d^)(l - + (1 - d^)r®' + (2 - d^ - d^ + 2d^°)r^ + 

' (1 - d^)fl + dS)r^ . 2d^(l + d^)r 2 - d^fl - d^) 


^27t 




Ig = / X cos'^G 008^^50 d 0 


4(l.d2)(r^.d2)(r^-d^)(r^-l)r^ U 
2 d^(l + d^)r^^ + (1 + d^°)r® + 2(1 - d^^)r^ + 

(1 - d^)(l - d^°)r^ - 2d^(l - d®)r^ - d^(l + d^°) 


2(1 + d2 + 2(1 - dMrl^ - 
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p2n 

= J X cos 0 cos^5® 
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Il6 


liy 


Il8 = 
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p 2n 
'^0 


X sin 0 cos 36 cos ^9 dd 


(: 


1 + 


^(1 + d^)(r^ + d^)(r^ - d^)(r^ - l)r 

2(1 - d^)r^2 + (1 + + (l - d®)r^ - (2 + d^ + d^ + 

S + (1 - d'*)(l * dS)!** + 2d‘‘(l + d«)r2 - d‘*(l - dS)] 


2d^0jr 


p 2jt 

■ / X sin 9 sin 3^ cos^0 d0 
^0 


jtr 


1 + d^)r^° + (1 


4(1 . d 2 )(r“ . d=)(r'* - . 1) 

(1 t d^)(l + d'*)r® - (1 - d®)r'* + d'*(l + d^)r^ + d'*(l 


d“)r8 


- d" 


n 2it 

^0 


X 'sin e sin 30 cos^30 d0 


1 - d^)r^^ - + d^]r^® +' 


l|-|l + d^)[r^ + - d^)(r^ - ijr^ 1— 

(l * d‘*)(l . d«)r« - (l - dl2)r'* - d“(l . d«)r^ * d'*(l - Jj 
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I-jQ = / X sin 0 sin 3^ cos 5^ 
Jq 


kil + d^)(r^ + d^)(r^ - d^ 


)(r^ - l)r^ 


2 1 - 




8 

d ll - d /r + 


1 - dN 1 - d^° 


\ 6 / ' , 16 \ 4 ,4(t ,10) 2 

j.T -. yl - d Jr - d \1 + d Jr + 


dMi - d"2 


I 20 = / 4 sin 9 sin 3® cos 9 cos 3© <3.0 

J 0 


4(l + d^)(r^ + d^jfr^ - d^ 


l(r'» - l) 


1 + d^jr® 


1 + d^Ul + d^Jr^ + d^fl + d^ 


1^^ = J \ sin 9 sin 3® cos 0 cos 50 d0 


41l + d^lfr^ + d^j (r^ - d^l/r^ - ijr 


, ,2] 14 

1 +, d IT 


1 + d^)(l + d^)r^° - (l - d®)r® + d^(l + d^jr^ + 


1 - d^)(l + dM^r^ - d^(l - d®. 
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X sin 0 3 ® cos 30 cos ^9 6.9 


4(1 + d^] (r^ + d^)(r^ - dM(r^ - 


1 + d2)rl6 - (1 - d^). 


<i^(l + d^)r^2 ■'■ (^ “ d^^)r^ - d^(l - d^)r^ + 


dMl ^ d^^ 


I 22 - J ^ sin e sin cos^G d© 


4(1 + d2)(r^ + d^)(r^ - d^) (r^ - l) 


1 - d^Jr^® + (1 + d^Jr® - 


1 - d^lr^ - 1 + d^)(l + d^jr^ + d^(l - d^jr^ + d^(l + d^ 


*-pk = / ^ sin G sin 5© cos 30 d© 

^0 


4(1 + d^Vr^ +'d.^Vr^ - dM(r^ - l] 


1 .+ d2)rl6 . (1 . 


1 + 6 ^) {2 - d^ + 2d^)r^^ + 2(1 - d®)r^° + 2d^(l + d^)r® - 


1 - d^)(l + 3d^ + d®)r^ + (1 + d^^)/ + 


d^(l - d®)r^ - d^(l + d^°' 
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‘■25 


p 2it 

J 0 


\ sin 6 sin 5® cos 5© 


# 

+ d^) + d^^l - d^jr^® + ^1 - d^^jp^ - |^1 + d^8)r*^ - 

d**(l - dl2)r2 + d^l + dl**) 


n2jt 

^26 “ / ^ ^ ^ 

Jq 


4(^1 + d^)(r^ + d^)(r^ - d^)(r^ - ijr^ l( 

1 + <3-^)[l + d^)r^*^ + ^1 - d®)r® + d^(l + d^)r^ - 


(1 - d^)(l + d^)^r^ + d^(l - d®) 


p2n 

^27 ~ ^ 50 cos 9 cos 50 d0 

^0 


4(1 + d2)(r'^ + d2)(r^ - d^)(r'^ - l)r"^ 

+ d^)^l + d^°)r^ + d^(l + d^°) 


(1 + d^°) 


r'" - 



32 


MCA TN 2326 





MCA TN 2326 


•33 


= J \ sin 30 sin ^0 cos^0 d0 


4(1 + d2)(r^ + d2)(r^ - d^) (r^ - ijr^ I- 


(l + d^)r^^ + - 2(1 - d^)rl2 


1 + d^)(l -,d^) - (1 - d®)r® - (2 + d^ + d^ + 2d^°)r^ 


1 - d-" Kl + d®)r^ + 2d^(l + d^)r^ + d^(l - d®) 


[32 = J ^ sin^3^ cos Q cos 3^ ^0 


k[l . . d^)(r^ - l)r^ " 

(l + d^)(l + d^)r^ - (l - dl^jr^ + d^(l + d^)r2 + d^(l - d^ 


1^2 = / 4 sin 5® cos 9 cos 3® d^ 


2 oA ^4 Vi6 

4(l.d2)(r^.d2)(rV.,4j(,4.ij,6.|^^ 

2(1 + d^)rl^ - (1 - d8)rl2 . (1 ■+ 2 d^ + 2 d^ + ^ 

d^(l - d^)r8 . (1 - d^)(l - dlo)r6 . (1 . + 


d^(l + d^'^jr^ + d^(l - d^^ 
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p 2jt 

= J X sin 30 sin . 50 cos , 6 cos 30 d0 
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X sin^0 sin^5® d® 


4^1 + 


- dMfr^ - l]r^ 


2(1 + d^)r^^ 


2(1 - dMr^^ - 2d^(l + d^)r^2 + ( 1 + d^°j 


j r® - 2 (1 - d^^j r^ + 


1 - d^lfl - d^^Jr^ + 2d^(l - d®)r^ - d^ (l + d^^J 


X sin 9 sin-^30 d0 


4(1 + d^ )(r^ + d^lfr^ - d^jfr^ - ijr 


3(1 - d^Jr^^ 


3(1 + d^)r^® • (^ ■ • 3<i^ + d^°)r^ * ~ + 


. j10i^6 


d^fl + d^lr^ - d^fl - d® 


= J X sin 0 sin 3^ sin ^6 dG 


4(1 + d2 + d2 '- d^ - llr 


1 + d^jr^^ - (l - d^jr^^ + 


1 + d^)(2 - 3d^ + 2d^ _ 2(1 - d®)r^° - 2d^(l + d^jr® 


1 - d^)(l - d^ + d®)r^ + (1 + d^Mr^ + d^fl - d^Jr^ - d^fl + d^^ 
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^53 ~ I ^ cos 50 

^ 0 


' ^ ^ (. hM^20 

4(1 . d2)(r“ * d^jfr** - d'*)(r‘* - l)r“ Li" ' ' 

(1 + al0)rl'* * (1 + d'*)(l + dl0)rl0 - d'*(l + dlo)r6 

I, , 2 oV 4 , 4 /, .leTI 


d Jr + d 1 - d 


\ sin.^50 cos ^30 d0 




2(1 + d^); 


1 + d^)(l + d^ +'d^)r^^ + 2(1 - d®)r^^ + d^(l + d^)r^^ 


2 1 - d^ 1 + d 


i^) r^° + 2d^(l - d®)r^ - (1 + d^®)r^ + d^(l + d^M 


= J X sin 5^ cos 30 cos ^9 d0 


4(1 + d2)(r^ + d2)(r^ - d^)(r^ - ijr^O 


1 - d^W20 + (1 + . 


3 )r 6 - (1 . d20)r^ + d"^(l - dl6)] 



NACA TN 2326 


4l 






42 


MCA TN 2326 




MCA TN 2326 


^62 ~ ^ sin^30 d.0 

'-'0 


3(1 + 


4 (1 + d 2 )(r^ + d 2 )(r^ - d^) (r^ - ijr"^ Ll^ ^ ~ ‘ 

3d2.^1 + d2)r^2 _ _ d^)r^ + -^1 + d^^) r^ 

d^(l + d^°) 


153=/' \ sln^30 sin 50 d0 

Jq \ 


(3 - d^ - d^ + 3dl0)rl0 + 3d^(l + d^)r6 + (l - d^^jr^ - d^(l - d^^) 

p2n 

/ X s±n^36 sln^^e d0 


21+ d^]r^° + 


4(l + d2)(r^ + d2)(r^- d^)(r^-l)r8 ' 

(1 + - 3d^ + d^)r^8 - 2^1 - d^jr^^ - d^ ^1 + d^jr^^ _ 

2(1 - d^)(l + d®)r^° + 2d^(l - d®)r^ + (1 + d^^)r^ - d^(l + d^^) 
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APPENDIX B 


THE Yn EXPRESSIONS 


After integration the functions 


in terms of tan 


-1 .-1 


and log. 


1 - d' 


Yn 

2 


can be put into closed forms 
where loge indicates natural 


1 + d^ 

logarithms. However , for small values of d, the computation 
eventually involves the differences of values of the same order of mag 
nitude. It is therefore necessary to expand these functions into infi 
nite series in terms of d. In such forms, the functions Yn are as 
follows : 


Yl = - ■ 2 2 + 1.333333333 - 0.200000000d^ + 

0.3l^2857ll^•2d^ - 0.053968253d^ + 0.202020202d® - 0 . 024975025d^® + 
0 . 1435897 ^ 3 ^^^ - O.Ol4379085d^'^ + 0 . Ill455107dl^ + . . . 


Y3 + loge 2 + 1.000000000 - 0 . 500000000d^ + 

4 (1 + di) (1 - d^ ) ' 

0.833333333d^ - 0.083333333d^ + 0 . 450000000d® - 0.033333333d^° h 


0.309523809d^^ - 0 . 0l7857l42d^^ + o. 236 iiiiiid^^ + . . . 


■Y5 = - —j—- - h logg 2 + 1.000000000 - 0.333333333d^ + 

2(1 + d ) (1- d j \ ^ 

0.533333333d^ r- 0. 00952 3809d^ + 0.253968253d® - 0.002020202d^° + 

0.l67832l67d^^ - 0.0007326011^^ + 0.1254901961^® + . . .] 


/ 
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Y7 = — ; -737 ot(-D - loge 2 + - 1 . 500000000d^ - 0 . l6666666Td^ 

‘ 4(1 + d 2 ) 2 (i . d 2 )\ 

0.583333333d® - 0.050000000d^° - 0.366666667d^^ - 0.023809525d^^ - 
0 . 267857 l 42 d ^^ + . . .) 


Yn = 


2(1 + d^) t- - 


(- i + I loBe 2 + d 2 . 1 . 333333333 d^ - 


0 . 466666667d® - 0.342857l42d® - 0.079365080d^^ - 0 . 202020202d^^ 
0.032l67833d^^ - 0.l43589743d^^ + . . 


Yii = E - 


4(1 + d^) (1 - d^) 


+ logg 2 + d^ - 1.500000000d^ + 


0.833333333d® + 0.4l6666667d^° + 0.450000000d^^+ 
0.l33333333d^^ + 0.309523809d^® + . . .) 


•-1 


D + log 2 + 1.500000000 - 0 . 333333333 d + 


4(1 + - d^) 

0.583333333d^ - 0.ll6666667d® + 0 . 366666667d® - 0.059523809d^® + 


0.267857l42d^^ - 0.036lllllld^^ + 0 . 2lllllliid^® + . . .) 


•-3 


2(1 + d^) (1 - d^) 


I - I logg 2 + 1.533333333 -_0.l42857l42d^ + 


o.253968253d^ - 0.0591630601^ + 0 . 1678321671^ - 0.0324786321^0 


0.1254901961^^ - 0.0205460161^^ + 0.1002506251^® 


+ . 
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Y_5 =-77- ttIt + loge 2 + 1.833333333 - 0.250000000d^ + 

U(l + d2)^ (l - d2) \ 

0.ii-50000000d^ - 0.ll6666667d^ + 0.309523809d® - 0.067857l42d^® + 
0.236lllllld^^ _ 0.044444444d^*^ + 0 . 190909091d^^ + . . .j 

Y_7 r-^7 r(| - I log 2 + I.676190475 - 0 . llllllllld^ + 

2(1 + d^) (1 - 

0.202020202d*^ - 0 . 056721056d^ + 0 . 143589743d® - 0 . 03458ll05dlO + 
0.lll455l07d^^ - 0.023327549d^^ + 0.09l097307d^^ + . . .) 


Y-9 r ^ rfr + logg 2 + 2.083333333 - 0.200000000d^ + 

.4(1 + d^j (l-d^)^ 

0.3666666671^ - 0.1095238091^ + 0.2678571421® - 0.0694444441^° + 
•0.2111111111^^ - 0.0480519491^^ + 0.1742424241^® + . . .) 

^-11 = - — — It -r(| - I log 2 + 1,787301587 - 0.0909090911^ + 

2(1 + i2) (1 _ d^j ' ® 

0.l67832l67d^ - 0.05268065.2d® + 0.125490196d® - 0.034532031d^° + 
0.l00250625d^^ - 0.02443064od^^ + 0.083478260d^® + . . ,j 

^-13 “ ~ ~ T27 ~\r ^ 2*2®3333333 - 0.l66666667d^ + 

4(1 + d2) (1 - d^) ® 

0.3095238091^ -0.1011904751® + 0.2361111111® - 0.0682539691^° + 
0.1909090911^^ - 0.0492424241^^.+ 0.1602564091^® + . . .) 
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Y-15 = - -r^T ^(1 - I log 2 + 1.878210678 - 0.076923076d^ + 

2 (1 + d^j (1 - d^j 

0.l43589743d^ - 0.048567120d^ + 0.11l455i07d® - 0.033583958d^° + 
0.091097307d^^ - 0.024654731d^^ + 0.077037037d^^ + . . }j . 

Y_17 r|- -fr +-logg 2 + 2.450000000 - 0.l42857l42d^ + 

4(1 + d^) (1 - d^j 

0.267857l42d'^ - 0.093253969d^ + 0.2lliiiiiid® - 0.065909091d^° + 

/ 

, 0.174242424d^^ - 0.049l45298d^^ + 0,l4835l647d^^ + . . !) 

0.i25490i96d^ - o.o4478844od^ + 0.100250625d® - 0.032273778d^*^ + 
O.08347826od^^ - 0.024405459d^^ + 0.071519795d^^ + . . .j 

Y_21 = - — r|- r(l> + logg 2 + 2.592857142 - 0.125000000d^ + 

4(1 + d2r(l - d2) ^ 

0.236lllllld^ - 0.086lllllld^ + 0 . 190909091d® - 0.0631313l3d^*^ + 
0.l60256409d^^ - 0.04835l647d^^ + 0.138095238d^^ + . . .) 

Y-23 = - — r(| - I log 2 + 2.021800421 - 0.0588235291^ + 

2(1 + d^j (1 - d^)' 

0.1114551071^ - o'.o4l427096d^ + 0.09l097307d® - 0.030846682d^° + 
0.0770370371^^ - 0.023900748d^^ + 0.06674o822d^^ + . . .) 


•-19 


2(1 + d^) 


d - I logg 2 + 1.955133754 - 0.066666667d^ + 
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Y = - — -i- -(d + log^ 2 + 2.717857142 - O.llllllllld^ + 

4(1 + d2)^(l - d2)^ 

0.21111lllld^ - 0.079T9T980d^ + 0 . 174242424d® - 0.060256409d^° + 
0.l4835l647d^^ - o.o47l86l47d^^ + 0 .l 29 i 66667 d^^ + . . .) 


2(1 + d2)^(l - d^) 


I - i log„ 2 +2.080623951 - 0.0526315781^ + 


0.1002506251^ -.0.0384657291^ + 0.08347826018 - 0.0294179901^° 


0.0715197951^^ - 0.0232625621^^ + 0.0625610941^8 + . , . 

where 



E = (logg) 

r=oo 

For 1 ^ = 0.075, the values of Yn are as follows: 

Yii = -0.163282840 - D' + E 
Y9 = 0.193584561 - D' 

, Y7 = -0.146562469 - D' 

Y5 = -0.295349926 - D' 

Y3 = -0.388303682 - D' 

Yi = -0.455432461 - D' 

Y_i = -0.507831136 - D' 

y_2 = -0.550751149 - D' 

Y_^ = -0.587074905 - D' 

Y_rj = -0.618549324 - D' 
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where D ' 


Y_5 = -0.646311118 - D' 

Y.ii = -0.671140539 - O' 
y_i3 = -0.693595951 - 0’ 

Y_^^ = -0.714090490 - D' 
Y-17 = -0.732938314 - D’ 
Y.19 = -0.750383821 - D' 

Y_2^ = -0.766620856 - D' 
Y-23 = -0.781805725 - O' 
Y_2^ = -0.796066296 - D’ 
Y_^rj = -0.809508489 - D' 

= 0. 233873678 [loggd - 
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APPENDIX C 

THE EXPRESSIONS 

The Kji expressions encountered in the calculation of flow past a 
curved surface by the variational method are as follows: 

= 0.059999539U^ + 0.l8295l490UAii^ - 0.28902V707UA11A13 - 0.099215517UA]^lAj^5 + 
0.0395^^924l^UAllA3l - 0.055l*^59*tOUAllA33 + 0.01l6870MtUAj^lA5j^ + 1.52l97Mt07UAi32 - 
1.708453568uA2^3Ai5 - 0.039032242UA3^3A33^ + O.462318173UA13A33 - O.OI0547763UA13A51 + 
1|..062796976 uAi 5^ - O.OOO252462UA15A31 - 0.203909l65UA3^5A33 + 0.000978600UA]^5A52. - 
0.0003874 o 6 uA 3 i^ •- 0.015012518 uA 32^A33 + 0.003802360UA33^A5 i + 0.03705372IUA332 - 
O.OO8I34189UA33A51 - 0.000749111UA5i^ + 0.033220706Aj^ 1^ + 0.022268910Ai3_^Ai3 + 
0.016688092 Aii %15 + 0 - 0.00ifl91746Aii2A23 + 0.010268355An^A5 ]^ + 
0.586559872A11A132 - 0.052992785A11A13A15 + 0.030572932A11A13A31 + 
O.196566992A12.A13A33 + 0.0178i6638A2^iAi3A5i + 1.590352068AiiA3^5^ + 

o 

0.027583320A11A15A31 + 0,026348o6iAiiAi 5A33 + 0.0l67973WfAnAi5A5i + 
0.007609784 AiiA 3 i^ + O.OO 8118406A11A31A33 + O.OO77219WA11A31A51 + 

0.01940l^931AllA33^ + 0.00503720 LAjj^A33A5i + 0.00214761AA]jlA5i^ + 0.3553^^8460A]^33 + 
1-998793077 Ai 3 ^Ai 5 + 0.126d29525Ai3^A3i + 0.1927ll803Ai3^A33 + 0.05^^55^*80^3^ 3^A5i + 
2,546743570 Ai 3A3^5^ + 0.093686l39Ai3A]^5A3i + O.725826526A13A15A33 + 
O.066165622A13A15A51 + 0.007521368 Ai 3A3 i^ + 0.0496 o 6 oi 4 Aj^ 3A3 j^A 33 + 
0.008865088A3^3A3 iA 53^ + 0.038026476A13A332 + 0.021958499Ai3A33A53^ + 

0.002747290A]_3A5 i^ + 1 .i 82535859 Ai 5^ + 0 . 35^628l07Ai5^A3i + 0.493373000Aj^5^A33 + 
0.158881*518 Ai 52 a ^1 + 0.006593000A15A312 + o.026235560Ai5A33^A33 +’ 

0.008096717A3_5A3iA5i + 0.068307146Ai5A33^ + 0.017553086Ai5A33A5j^ + 

O.OO2597463A15A5X2 + 0.001063060A3 x 3 + 0.002385806A3 x^A 33 + O.OOl8l66lOA3x^A5x + 
0.005491591A3xA33^ + 0.003005606A3XA33A5X + 0.0011o6980A3xA5x^ + 0.002670255A333 + 
0.0028 i 6774A33^A5 x + 0.000987380A33A5x^ + 0.000237098A5 x^ 
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Kg = - 0.130908999 u3 - 0.1kk^l233^f^^^ + 3.o439W8l4UA^j^A^^ - 1.708453568UA^^A^^ - 
0.0390322l^2UAllA3l + O.462318173UA11A33 - O.OIO547763UA11A51 + 3 . 21377H8UAi3^ + 
11.8579^352UAi3Ai 5 + 0.346^76o6irtJAi3A3i + I.03127970I1UA13A33 + 

0.0917699'^6uAi3A5i + 7.73l273421UAi5^ - 0.3^98o8o85UAi5A3i + I.726043683UA15A33 - 
O.146059989UA15A51 + 0.000692852UA3i^ + 0.058004535UA3j^A 33 + 0.002368092UA3iA5j^ + 

0. 0872000 38 uA 33^ + 0.060930600UA33A51 + 0,001909217UA5 j^^ + 0.007422970Ai]^^ + 

0.586559872A^j^^Aj^ 2 - 0.026496392A^j^^A^5 + 0.015286'466 Aj^^^A2j^ + 

0 .09Q2Q3hsep^^k^^ + 0.008908319Alj^^A^^ + 1.066045380A^Aj^^^ + 
3.997586154 A]^iAi3Ai5 + 0.252059050A11A13A31 + o,3854236o6Ai3^Ai3A33 + 
O. 10910968OA11A13A51 + 2.564743570A11A152 ^ 0.093686139A11A15A31 + 
0.725826526a^j^A^^A^^ + 0.066i65622A^^A^^A^^ + 0.007521368a^j^A2^^ + ; 

o.o496o6oi4Aj^^A23^a^2 + 0.008865088a^3^A3^A^^ + 0.038026476Aj^ 3^A23^ + 
0.021958499A]^3^A22A53^ + 0.p02747290A^^A^^^ + 3.60928926 oa^ 2^ + 8.7l9639l40A^3\5 + 
0.208248666 a^2^a^^ + i.8i174i922A^2^A22 + o.o85998486a^^^a^j^ + 

25.4i4218o49Ai3Ai 5^ + 0,920228l70A3^3A]^5A33^ - 3.000980028A3^3A3^5A33 + 
0.455104944 a^^A^^^^ + 0 .03orjk6x&\k^^^^ + 0.086959326a^^A23^A2^ + 
0.042256i62Ai3A3j^A5i + 0.324078519Ai3A33^ + 0.040162146A13A33A51 + 
0'.012506092Ai3A5]^^ + 10.87146o415A3^53 + o.504585238A]^5^A3j^ + 4.344o47335A3^52;^33 + 
0,210736557A3^5^a^3^ + 0.034424058a^^A2^^ + 0.168542116a^^A23^A23 + 
0.04555146ia^5A3j^A^j^ + 0 .2Q9&9hQ36k^^k^^ + 0.040344011A^^A33A^^ + 

0 .Ql%k'^S^lky^P^-^ + O.OOIO33455A313 + o.008o85122A33^2^33 + 0.001927351A3]^2a5i + 

0.009687255A3j^A332 + o.009037828A3^A33A^^ + 0.001182158A3j^A^^2 + 0.020369802A333 + 

0.004917098A33^A^^ + 0.002764734A33A^j^^ + 0.000247234a^^3 



0.009387226U^ + 0.0197T4622UAii^ - 0.0390322l+2UAiiAi3 - O.OOO252U62UA11A15 - 
0.0p07748l2UAj^3^A2j^ - 0.015012518UA^^A23 + 0.003802360UA^j^A^^ + 0. 1732380 32UAj^^^ - 
0.3'f98o8o85UAi3Ai^5 + O.OOI385705UA13A31 + 0.05800I+535UA13A33 + 0.002368092UAi3A5j^ + 

0A2787386J(UAj^^^ + 0.020248502UA^^A2^ - 0. 0359^720 3UA^^A^2 + 0.007519578uAj^^a^j^ - 
0.002531840UA2j^^ - 0.0'04441)-278UA2^A22 - 0. 000^6698^1 ! + 0.0037l^9007UA22^ - 
0.0020907U9UA23A^3^ - 0.00107891 hua^^^ + 0.007876620A^^^ + 0.015286J*66Aj^^^Aj^2 + 

0.013791660Aj^^^;^^ + 0.007609784a^^^A 2^ + 0. 00405920 3Aj^^^A^2 + 0.003860974 Aj^j^^A^j^ + 

0.126029525Aj^j^A^2^ + 0.093686139A^^A^^A^^ + 0,015042736Aj^j^A^2A2^ + 

0. 04960601 4a^j^Aj^2^33 0-00®®^5088a^^A^2A5i + 0.354628107A^^j^a^^^ + 

0.013186000Aj^j^A^^A^j^ + 0.026235560A^^A^^A^2 + 0.OO8O96717Aj^3^Aj^5A^j^ + 

0.003189180Aj^j^A^^^ + 0.004771612A^3^A23^A^2 0.003633220Aj^j^Ajj^A^^ + 

0.005491591A^j^A22^ + 0.003005606A^^A2^A^j^ + 0.001106980Aj^^A^j^^ + 0.0694l6222A^^^ + 
0.460114085Aj^2\5 + 0.039746i 67A3^3^A2^ + 0.04347-9663Aj^^^A23 + 0 . 02 ll 28 o 8 lAj^ 2 %i + 
0.504585238a^ A 2 + 0.068848115A A A + 0.i685421i6a a A '+ 

13 15 13 15 31 13 15.33 

0.04555146ia^2\^A^^ + 0;003100365A^3A2^^ + 0.016170245A^^A^^A23 + 
O.OO30547O3A^2A3iA5i + 0.009687255A^3A^2^ + O . 009037828a^^a^^A^^ + . 
0.001182158a^2A^^^ + 0.230447196a^^^ + 0.120600559A^^^A2j^ - 0.260365077A^^^A^3 + 

0.06734O386a^^^A^^ + 0.002708871A^^A3^^ + 0.01538l835Aj^^A^^A32 + 

. 2 

O.OO329776IA15A31A51 + O.017641290A15A33 + 0.01007441 5A]^^A33A5i + 
0.001024491A^^A^^^ + 0. 000551 450A3^^ + 0.001323410A3^^A33 + 0.001028092A^^^A^^ + 

0.002235226A3^A3^^ + O.OO170646lA2^A33A^^ + 0.000672957A3 j^A^j^^ + 0.000845763A^2^ + 
0.00134i68iA33^A5i + 0.000569815A33A^]^^ + 0.000153120A5jl^ 
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3 2 

K4 = -O.OI6023195U - O.02757297OUA11 + O.462318173UA11A13 - 0,203909165 UAj^iAi 5 - 

0.015012518qa^j^A2j^ + 0.074i 07442UA^A23 - 0.008l34l89UA^j^A^j^ + 0.5l57l3852UAj^^^ + 
1.726043683UAj^2Ai 5 0.058004535UAj^2A2j^ + 0.174400076UA^jA22 + O.06O93060OUAj^2A5i + 

, 1.217327709UAi 5^ - 0.0359^^7203UAl5A3l + 0.23575*H37UAi5A33 - 0.0212206JnUAi5A5i - 

0.002222139UA3i^ + 0.007‘^98011(UA3lA33 - 0.0020907^^9UA3lA5l + O.Ol41H9153UA33^ + 
O.OOI269098UA33A51 - 0.00050938toA5i^ - 0.0013972i^9All^ + 0.098283496 Aii%3 + 
0.01317*tO31Ajj^^Aj^^ + 0.004059203 A^j^^A2j^ + 0.019'K)4931A^^A22 + 0.002518601A^j^^A^^ + 
0.1927U803A^^A^2^ + 0,725826526a^j^A^ 2A^^ + O.Ol^96o6oll^AJ^J^A^2A2J^ + 
0.076052952 Ajj^A^ 2A33 + 0.0219581t99A^3^A^3A^j^ + 0.493373000 A^Aj^^^ + 
O. 02623556OA11A15A31 + O.136614292A11A15A33 + 0.017553086AiiA3^5A5j^ + 
0.002385806a^j^A2j^^ + 0.010983182A^j^A23^A^3 + 0.0030056o6Ajj^A^j^A^^ + 
0.008010765A3^iA33^ + 0.005633548Ai3^A33A52^ + O.00O987380A]jl^1^ + 0.603793974Aj^3^ - 
1,500490014Aj^2^Aj^5 + 0.o 43479663A^2^A^^ + 0.324078519Aj^2^^33 + 0.02008l073A|^3^A^3^ + 
4.344 o47355A3^3A^5^ + 0.i68542116a^2Ai5A3i + 0 + 

0.040344011A^2Aj^^A^^ + 0.008085122A^2A3i^ + 0.019374510A^2A3iA33 + 
0.009037828 Aj^2^2iA5i + 0,06i109406A3^2A33^ + 0.009834196Aj^2A33A5i + 
0.002764734A3^3A^3^^ + 1.932393855A^^3 - 0 .2(X)Z6yy9 ^^ + 0.8o6221543A^^^A23 + 
0 .o6o527562A^^^A^3^ + 0.007690918A|^5A2j^^ * O.oyys&Z^QOk^^k^^k^^ + 
0.010074415A3^5A3j^A53^ + 0.056866015Aj^^A 33^ + O.Ol8945942Aj^5A33A5j^ + 
0.003407649A^^3^^ + 0.00044 i137A2j^^ + 0.002235226a^^^A 22 + 0.00085323lA23^^A^j^ + 
0.002537289A3j^A23^ + 0.002683362A2j^A 23A^^ + O.OOQ5698l5A23^A^j^^ + 0.004ll5250A23^ + 
0.001435584A22^A^3^ + 0.000864839A23A^j^^ + 0.000129993A^^^ 
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= 0.0092414310^ - 0.0496O7758uAi]^^ - 1.708453568UAnA2^3 + 8.125593952UAj^lAi5 - 

0.000252462UA^^A2^ - 0.203909165UA^A^2 + 0 . 000978600UA^^A^^ + 5.9286l99l8UA^2^ + 

15.462546842UAj^ 2^^ - 0.3498080850/^2^31 + l.726o436830Aj^2A33 - 

O.146059989OA13A51 + 10.786i660480Ai5^ + O.855747728OA15A31 + 2.434655418OA15A33 + 

0.1908i89700A^^A^^ + 0.0101242510A2i^ - 0.0359472030A^^A23 + 0.0075195780A2iA^l + 

0.1178T7054oA23^ - 0.0212206410A23A^1 + 0.0031449810A^^^ + 0.005562697A^^^ - 

0.026496392A^i\3 * 1.590352068a^^^A^^ + 0.01379l660A^^^A2j^ + 0.013l74031Aj^l^A22 + 

0.008398672A^1^A^1 + 1.998793077AnAj^3^ + 5.l29487l40A^lAi2Ai5 + 

0 . 093686139^1 T A-) + 0.725826526Aj^ 2,A^2A^^ + 0,G66l65622AQ^2^Aj^y^2^ + 

2 

3.5^7607577A2iAj^^ + 0.709256214Aj^ 2.-A-]^5A21 + 0,9867^000A^Aj^^^2 + 

0.317769036a^iA^iA^ 1 + 0.006593000/^j^A2i^ + 0.026235560A3j^A2iA23 + 

2 

0.008096717A11A31A51 + 0.068307146A11A33 + 0.017553086A11A33A51 + 

0.002597463Aj^^A^^^ + 2.906546380Aj^3^ + 25.414218049/^3^^ + 0.46 oii4o85A3^2%i ‘ 

1.500490014Ai3^A33 + 0.227552472Ai3^A5i + 32 . 6l438l245Ai3Ai5^ + 

1.009170476a^2Ai5A2i + 8.688094670/^^2^15^33 + 0.42l4731l4/l|^2Ai5A5i + 

0.034424058a^^A^^2 ^ o.i68542116a^^A^^A^^ + 0.04555l46u^^A^^A^^ + 

O. 289894836 A A + 0.040344011A, A A , + O.OI 564565 OA A 2+ 

13 33 51 13 51 ' 

36.413233786/1^^3 0.691341588Aj^5^A2i + 5.797l8l565A3^5^A23 + 0.287053533A^5^A^1 + 

O.12060O559Ai5A33^^ - 0.520730154/1^^21^3 + 0.134680772Aj^5A23^a^j^ + 

0.8o6221543Aj^^A 22^ + 0.121055124/1^^^23*51-+ 0.04l784273A^^A^^^ + 0.000902957A2^3 + , 
0.007690918A2i^A 22 + O.OOl64888lA2i^A^3^ + 0.01764l290A2iA22^ + 

O.OIOO74415A31A33A51 + 0.001024491A3iA5i^ + 0.018955338A33^ + 0.00947297U33^A51 + 
0.003407649A23A^1^ + 0.0002l4347A^l3 
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^ 2 

Kg = O.OO3281673U + O.OO5843522UA11 - 0.0105J+7763 UAiiAi3 + O.OOO97860OUA11A15 + 

O.OO380236OUA11A31 - O.OO8134189UA11A33 - O.OOI498222UA11A51 + 0 ; 045884973 UAi 32 - 
0 . 046059989 UAj^ 3 Aj ^5 + 0 , 002368092 UA^ 2 ^^ + 0 . 060930600 UA^jA 23 * 9 . 0038 l 8433 UA^ 2 ^j^ + 

0 . 095409485 UA^^^ + 0 .001^19^ieVA^^A^^ - 0 .0212Z06kVJA^^A^^ + 0.00628996lUA^^A^^ - 

2 '■ 2 ' 

0.000283494uA32^ - 0.002090749UA31A33 - '0.002157882UA3iA52^ + O.OOO634549UA33 - 

O.OOi6i876tuA2^A^3^ - 0.000752348ua^^^ + 0 . 003422785A^ + 0.008908319A^3^^a^3 + 

0.008398672A^^A^^ + 0.003860974a^]^^A^^, + 0.0025l860lA^^A^2 + 0.002 i47644a^^^a^^ + 

0.054554840A A ^ + 0.066165622A A A + 0.008865088A A A + 

11 13 - 11 13 15 11 13 31 

0 . 021958499 A^A ^^23 + 0 . 005494580 A^^A^^A^^ + 0 . 158884518 a^^a^.^^ + 

0.008096717A;l^A^5A2j^ + 0.017553086a^3^A^^A22 + 0.005194926a^A^^A^^ + 
0.0018 i 6610A^A2^^ + 0.00300^6o6A^^A^^A^^ + 0.002213960A^A2j^A^j^ + 
0 . 0028 i 6774 A^A^ 2 ^ + 0 . 001974760 A^j^A 22 A^^ + 0 . 0007 U 294 a^A^j^^ + 0.028666l62A^2^ + 
0 . 227552472 A^ 2 ^A ^5 + 0.021128o8lA^2%i + 0 . 02008 l 073 A^g %3 + 0 . 012506092 Aj^ 2 %i + 
0 . 210736557 A 3 ^ 3 A^ 5 ^ + 0 .o 4555 l 46 lAj^ 2 Ai 5 A 3 i + 0 . 04034401 LA^ 2 A ^^23 + 
0.031291301A^2^^A^^ + 0.001927352A^3A2;l^ + 0.009037828A^2A23^A22 + 
0.002364316A^3A23^A^j^ + 0.0049170y8A^2A33^ + 0 . 005^29 468a^2^2 3^51 + 
0.000741701A^jA^^^ + 0.0956845HAj^5-^ + 0.067340386A3^5^A2j^ + 0.060527562Aj^5^A23 + 
0.041784273Aj^5^A^j^ + O.OOl64888lA^^A23_^ + 0.0100744l5A^5A2j^A22 + 

0.002048982A^^a^^a^^ + o.oo9472971A^5A22^ + o.oo6815299A^^a^2A^^ + 

2 3 2 2 
0.000643042Aj^5A5i + 0.000342697A31 + O.OOO853231A31 A33 + 0.000672957A33^ A51 + 

2 , 2 ‘ 

O.OO134168IA31A33 + 0.001139630A3]^A33A5i + 0.000459360A3JLA51 + 

Q ' 2 ^2 ^ 

0.000478528A33 + O.OOO864839A33 A^2. + 0.000389980A33A5]^ + 0.000108272A5]_ 
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TABLE I 


PARAMETERS Aij/ao FOR FLOW. PAST A CIRCULAR CYLINDER 
AT VARIOUS UNDISTURBED-STREAM MACH NUMBERS 

Q-«] 


M 

0 

A,, /a 
11/ 0 

0 



on . 

5i/“o 

A__ /a 
33/ 0 


^51/^0 

0.1 

0.001031 

.001035 

.001103 

.001103 

.001103 

.001105 

-0.0002557 

-.0002557 

-.0002541 

-.0002542 

-.0002541 

-0.0004301 

-.0004301 

-.0004301 

-.0004378 

-0.00001330 

-.00001331 

-.00001337 

0.000001180 

.000001180 

0.00001300 

0.2 

0.008605 

.008747 

.009347 

.009350 

.009346 

.009568 

-0.002199 

-.002198 

-.002140 

-.002141 

-.002135 

-0.003790 

-.003780 

-.003780 

-.004123 

- 0.0004862 

-.0004979 

-.0004981 

0.00004361 

.00004361 

0.0005734 

0.3 

0.03114 

.03278 

.03520 

.03518 

.03519 

.03540 

-0.008506 

-.008503 

-.007949 

-.007973 

-.007972 

-0 ’01514 
-.01491 
-.01489 
-.01803 

- 0.004733 

-.005045 

-.005051 

0.0004367 

.0004367 

•0.005229 

0.4 

0.08307 

.09430 

.1026 

.1018 

.1025 

.1038 

-0.02580 

-.02581 

-.02225 

-.02265 

-.02264 

- 0.04693 

-.04381 

-.04297 

-.06157 

-0.03156 

-.03722 

-.03727 

0.003155 

.003156 

* 

0.03036 

0.5 

0.1979 

.3242 

.3541 

.3431 

-0.09897 

-.09980 

-.06491 

-0.1657 

.05611 

-0.4292 

i 
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TABLE II 

MAXIMUM VELOCITIES q/U OVER A CIRCULAR CYLINDER COMPUTED BY 
VARIATIONAL AND RAYLEIGH- JANZEN METHODS AT VARIOUS 
UNDISTURBED-STREAM MACH NUMBERS 


\ 3. 
\u 

- 


Variational method 


' 

Raylei^-Janzen 
^method (third 
approximation ) 



y 

= 2 



7=1.5 

Mo\ 

One 

term 

Two 

terms 

Three . 
terms 

Four 

terms 

Five 

terms 

Six 

terms 

One 

term 

7 = 1.405 

7 = 2 

0.1 

2.0069 

2.0120 

2.0119 

2.0119 

2.0120 

2.0120 

2.0069 

2.0119 

2.0120 

.2 

2.0287 

2.05115 

2.0506 

2.0510 

2.0518 

2.0524 

2.0284 

2.0513 

2.0520 

.3 

2.0692 

2.1296 

2.1282 

2.1309 

2.1363. 

2.1364 

2:0681 

•2.1314 

2.1360 

.h 

2.1385 

. 2.2862 

2.2844 

2.2979 

2.3333 

2.3336 

2.1040 

2.2836 

2.3027 

.5 

?'2.2639 

^■2.8281 

^•2.8271 

^3.0754 



2.2494 

2.5707 

2.6312 
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TABLE III 


PARAMETERS FOR FLOW PAST A KAPLAN , BUMP AT VARIOUS 

UNDISTORBED-STREAM MACH NUMBERS 


M 

^11/^0 

^13 1^0 

^31/^0 

^33/^0 

% 5 /^o 

^l/^o 

0.5 

0.009459 

-0.003325 

-0.001799 

0.000558 

0.000180 

0.000156 

0.75 

0.05289 

-0.02119 

-0.01607 

0.01662 

0.02490 

-0.01327 

0.83 

0.10501 

- 0.04533 

-0.02624 

0.04201 

0.00799 

-0.07141 

0.90 

0.092073 

.20799 

.2554 

.2736 

-0.081943 

-.08626 

-.1064 

-0.1917 

-.2256 

0.1558 
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■Four parameters of the series taken. 
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TABLE V 

VALUES OF SIN 9 , SIN 30 , SIN 50 , AND 1 ^ 

I dz 

FOR VARIOUS VALUES OF X 


X 

sin 0 

sin 30 

sin 50 


—I 

dz 1 

0 

1.00 

-1.000 

1.000 

1.8500 

.1 

.99^ 

-.950 

.8628 

1.8434 

.2 

.9777 

-.8053 

.4909 

1.8226 

.3 

.9^92 

-.5732 

-.0300 

1.7953 

.h 

.9086 

-.2746 

-.5518 

1.7318 

.5 

.8544 

.0684 

-.9172 

1,6537 

.6 

.7850 

.4202 

-.9794 

1.5469 

.7 

.6963 

.7385 

-.6513 

1 , 

3997 

.8 

.5807 

.9589 

.0435 

i.1926 

• 9 , 

- .4186 

.9626 

.8326 

.87926 

.975 

.2119 

.6156 

.9609 

• 

45294 

1.00 

0 

0 

0 

0 
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TABLE VI 

PRESSURE COEFFICIENTS COMPUTED BY VARIATIONAL METHOD 
FOR VARIOUS UNDISTURBED-STREAM MACH NUMBERS AT 
VARIOUS VALUES. OF X ON BUMP 


\Cp 

X 

0 

II 

LTN 

• 

0 

II 

% = 0.75 

= 0.83 

0 

-0.168777 

-0.196976 

-0.278837 

-0.380572 

.1 

-.164025 

-.191428 

-.270078 

-.366493 

.2 

-.151114 

-.176112 

-.245300 

- . 327481 

.3 

-.130820 

-.151796 

-.206697 

-.266558 

.4 

-.101030 

-.117012 

-.148537 

-.187512 

• .5 

-.067709 

-.077888 ' 

-.095925 

-.102108 

.6 

-.030022 

-.033808 

-.033924 

-.017223 

.7 

.010174 

.012772 

.027712 

.058712 

.8 

.051714 

.060936 

.085342 

.118537 

.9 

.093315 

.107692 

.135733 

.155920 

.975 

.124466 

.143040 

.167266 

.161538 
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TABLE VII 

PRESSURE COEFFICIENTS COMPUTED BY VARIOUS METHODS 
FOR Mq = 0.83 AT VARIOUS VALUES OF X ON BUMP 


\Cp 

X N. 

Prandt 1 -Glauert 
method 

Von KanoEui-Tsien 
method 

Perturbation 

method 

Vaxiational 

method 

0 

-0.30245 

-0.324138 

-0.35564 

-0.380572 

.1 

-.294033 

-.314480 

-.34310 

-.366493 

.2 

-.270905 

-.288167 

-.30771 

-.327481 

.3 

-.234545 

-.247375 

-.25271 

-.266558 

.4 

-.181081 

-.188634 

-.18424 

-.187512 

.5 

-.121378 

-.124726 

-.10861 

-.102108 

.6 

-.053787 

-.054434 

-.03270 

-.017223 

.7 

.018234 

.018160 

.03827 

.058712 

.8 

.092715 

.090853 

.10222 

.118537 

.9 

.167303 

.161335 

.16275 

.155920 ■ 

• 975 

. 223151 

.212658 

.21232 

.161538 
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Angle, 0, deg 

Figure 2 . - Velocity distribution over a circular cylinder at Mq = 0.4 

computed by variational, Rayleigh-Janzen, and incompressible -flow 
methods. 




Maximum local Mach number 
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Figure 3 .- Maximum local Mach number at various undisturbed -stream Mach 
numbers computed by variational and Rayleigh-Janzen methods. 




Local Mach number at surface of cylinder 
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Angle, 9, deg 

Figure 4,- Local Mach number over surface of cylinder computed by 
variational and Rayleigh-Janzen methods for various undisturbed - 
stream Mach numbers. 




Method 

Variational 
Incompressible -flow 
Karm^n-Tsien 
Rayleigh-Janzen 





Pressure -coefficient ratio 
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Undisturbed -Stream Mach number, Mq 

Figure 6.- Pressure -coefficient ratio for circular cylinder -at various 
undisturbed -stream Mach numbers computed by various methods. 








Pressure coefficient 
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X 

Figure 8.- Pressure coefficients computed by variational and 

incompressible -flow methods for various imdistur bed -stream Mach 
numbers at various values of X on bump. 










